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Abstract 

In  recent  years  the  basic  structure  required  to  implement  the 
inverse  scattering  transform  in  1+1  and  2+1  dimensions  has  been 
clarified  and  extended.  Aspects  involved  with  fully  multidimensional 
problems  have  also  been  treated.  In  particular  the  inverse  scattering 
associated  with  various  multidimensional  operators  and  generalizations 
of  the  Sine-Gordon  and  self-dual  Yang-Mills  equations  have  been  studied. 
A  review  of  some  of  this  work  will  be  discussed  in  this  review. 
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^  The  Inverse  Scattering  Transform  (I.S.T.)  is  a  method  to  solve 
certain  nonlinear  evolution  equations.  There  has  been  wide  ranging 
interest  in  this  method  for  many  reasons^)  A  review  of  earlier  work  can 


be  found  in  [1],  A  surprisingly  large  number  of  physically  interesting 
nonlinear  equations  can  be  solved  via  1ST;  there  are  many  applications  in 
physics  including:  surface  waves,  internal  waves,  lattice  dynamics,  plasma 
physics,  nonlinear  optics,  particle  physics  and  relativity.  Mathematically 
speaking  the  field  is  also  quite  rich,  with  nontrivial  results  in  the 
areas  of  analysis,  group  theory,  algebra,  differqptial  and  algebraic 
geometry  being  used  by  various  researchers. 


differential  anq  aig« 
i .  FronOcffiP^poi  nt“  of 


vien  1ST 


allows  us  to  solve  the  Cauchy  problem  for  these  nonlinear  systems 

>-shaH  concentration  questions  in  infinite  space.  All  of  the  nonlinear 

A 

equations  discussed  below  arise  as  the  compatibility  condition  of  certain 
linear  equations,  one  of  which  is  identified  as  a  scattering  (direct  and 
inverse  scattering  is  required)  problem  and  the  other(s)  serves  to  fix 
the  "time  evolution1'  of  the  scattering  data.  <<r — 

In  one  spatial  dimension  the  prototype  problem  is  the  (KdV) 
equation 

u..  +  6uu  +•  u  *0. 
t  x  xxx 

The  KdV  equation  is  compatible  with 

xv 


vxx  ♦  U ( X  , t ) V 


vt  *  (Y+Ux)v  ■  (4*+2u)vx 


(1) 

(2) 

(3) 


i.e.  vxxt  *  vtxx  implies  (I).  Equation  (2)  is  the  time  independent 
Schrodinger  scattering  problem,  x  the  eigenvalue  (y  3  const,  in  (3)).  The 
solution  of  (I)  on  the  line:  -«<x<»  for  initial  values  u(x,ts0) 
vanishing  sufficiently  rapidly  at  infinity  is  obtained  by  studying  the 
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associated  direct  and  inverse  scattering  problem  of  (2)  and  using  (3) 
to  fix  the  time  evolution  of  the  scattering  data.  It  turns  out  that 
the  inverse  problem  amounts  to  solving  a  matrix  Riemann-Hi lbert 
boundary  value  problem  (RH8VP)  whose  jump  discontinuity  depends 
explicitly  on  the  scattering  data.  Calling  k»-k2,v(x,k)«u(x,k)e'’kx 
the  RH8VP  takes  the  following  form, 

(u4-uj(x,t,k)  *  u_(  x,t,t  (k))  V(x,t,k)  on  1 

wt-l,  lk|—  (4) 

where 

V(x,t,k)  .  r(k,t)  e2lkx,  a(k)  -  -k.  E»{k:kc«),  and  are  the 

limiting  boundary  values, as  Imk+Ot,  of  meromorphic  functions  in  the 
upper  (♦)  lower  (-)  half  plane.  (4)  may  be  converted  into  a  linear 
integral  equation  by  taking  a  minus  projection  and  the  potential  is 

reconstructed  via 

u(x.t)  .  -  1  ;  u(k.x.t.-k)  V(x,t,k)dk  (5) 

»  ax  j 

where  the  contour  is  taken  above  all  poles  of  r(k,t);  of  which  there 
is  at  most  a  finite  number,  k.  *  ir.t  «c  '0  j  *  1,***N.  The 

J  J 

scattering  data:  the  reflection  coefficient,  r(k.t)  evolves  simply 
in  time  2 

r{ k , t )  *  r(k,0)  e^’k  C  (6) 

The  above  scheme  may  be  extended  so  as  to  solve  a  surprisingly 
large  number  of  interesting  nonlinear  evolution  equations.  There  are 
two  scattering  problems  of  particular  interest  in  one  dimension: 
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(i)  Scalar  scattering  problems: 
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Via  an  appropriate  transformation  the  inverse  problem  associated  with 
(i),  (ii)  can  be  expressed  as  a  matrix  RH8VP  of  the  form  (4).  The 
potentials  Uj.q  can  be  shown  to  satisfy  nonlinear  evolution  equations 
by  appending  to  (i)  and  (ii),  suitable  linear  time  evolution  equations. 
One  then  finds  that  the  scattering  data  V(x,t,k)  evolves  simply  in 
time.  Well  known  solvabTe  nonlinear  equations  include  the  Boussinesq, 
modified  KdV,  sine-Gordon,  nonlinear  Schrodinger,  and  three  wave  in¬ 
teraction  equations.  The  reader  may  wish  to  consult  for  example  [2a-e] 
for  a  detailed  discussion  of  some  of  this  material. 

It  is  most  significant  that  these  concepts  can  be  generalized  to 
2  spatial  plus  one  time  dimension.  Here  the  prototype  equation  is 
the  Kadomtsev-Petviashvili  (K-P)  equation: 

2 

(u*  ♦  6uu  ♦  u  )  *  - 3o  u  (7) 

•  t  x  xxx  *  x  yy 

which  is  the  compatibil ity  equation  between  the  following  linear  prob¬ 
lems: 

av  ♦  v  ♦  u(x,y,t)v  *  0  (8) 


v.  ♦  4v  ♦  6uv  ♦  3(u  - o 

t  XXX  X  X 


Uydx ' )v  ♦  >v  =  0  (9) 


(y  *  const.).  We  shall  consider  the  question  of  solving  (7)  for 

2  2  2 

u(x,y,0)  decaying  sufficiently  rapidly  in  the  plane  r  =  x  ♦  y  ■*  »\ 

?  2 

Physically  speaking,  both  cases  o  *  -1  ( KP I )  o  s  «-l  (KPII)  are  of 

interest.  Whereas  KPI  can  be  related  to  a  RHBVP  of  a  certain  type 

(nonlocal;  see  ref.^)  KPII  turns  out  to  require  new  ideas.  Letting 

2 

/  ..  ikx  ♦  k  y/o 

v  «  ;i(x,y,k)e  7 


o  *Cp  ♦  i  Tj ,  -p  t  0.  Then  there  exist  functions  .  bounded  for  all 
x  ,y  satisfying  •  1  as  |k  |  »  <**.  However  such  a  function  turns  out 
to  be  nowhere  analytic  in  k,  rather  it  depends  nontrivially  on  both 
the  real  and  imaginery  parts  of  k(k»kR  ♦  ikj).  u  «  u(x ,y ,kR,kj ) . 

In  fact  u  satisfies  a  generalization  of  a  RHBVP  -  namely  a 
3  (08AR)  problem  where  u  satisfies. 
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“■*  *  u(*»y»Cg»^j)  v(*»y*kR*kj) 

V  has  the  struct 
iB(x#yfkRtkj « £q ) 


3t  2'3kfi  '  5kI 
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where  ±  .  i(-f-  ♦  i  J-)  and  V  has  the  structure 


Vfx.y.kfl.kj)  « 


sgn(kg)e 


2"I°rI 


I{kR.k,) 


»(«.y.kR.k,.«0)  ■  (»  *  2y  ^>(e0  -  k(l>  ■  -2(x  *  2y  4,lto 
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(11) 
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(10-11)  may  be  converted  into  a  linear  integral  equation  by  employing 
the  generalized  Cauchy  formula.  T(kR,kj )  is  viewed  as  the  "nonphysi¬ 
cal"  data,  (i.e.  inverse  scattering  data  or  inverse  data)  and  the 
potential  is  reconstructed  via 

(12) 


„(«,,)  •  ^  i  j  j  »(*.y.t0.kI)»(».y.k*.kI)«kR<ikI. 


The  basic  ideas  used  in  order  to  derive  these  equations  is 
as  follows.  We  convert  the  equation  for  u  *  u(x,y,k): 


tJUy  ♦  uxx  +  2ikux  -  u(x,y)u  *  0 


(13) 


into  an  integral  equation 

u(x,y,k)  »  1  ♦  G(u,u) 


(14) 


where 


G(f)  *  G*f  -  [[g(x-x‘ ,y-y ' ,k)  f(x' ,y' )dx' dy\ 

*  4 

the  Green's  function  kernel  being  givenby  (k*kR+ik, ) : 


(15) 


G(x ,y ,kp ,kj ) 


iUx+yy) 


dCdy 


(2 n)c  ( ion-£  -2k£) 

.  sg^(y)  fdCe’X5  +  C(p2k)y/o 


2*o 


3  (  -yoR(C^2CkQ))dC 


(16) 
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where  kg*kR  *  -  k[  and  J  (x)  s  '  lx>0»  Ox < 0 r 

R 

The  5  derivative  of  the  Green's  function  is  especially  simple. 


3G(  k  k  )  s  ii^)  iS(x,y,kR,kr 

-Ix.y.Kp.Kjj  e 


3/3C  =  |(ll^  *  1  and 

kT 

a(x,y,kR,kj)  *  -2(x+2yj— ))Cq. 

Taking  the  5  derivative  of  (14) 

^(x.y.kp.kj)  ’  jr!^f(x*x‘  V  »kR.kj)u(x'  ,y')u(x'  ,y'  .kp.k^dx'dy' 

♦  G(x-x' ,y-y' ,kR,kj)u(x' ,y' )^r-(xl .y' ,kR,kj )dx' dy '  (18) 

and  using  (17)  shows  that 

3u  sgn(kn) 

T(kRikx)  wfx.y.kp.kj)  (19) 

where  T( kR , k j )  *  e"’8^x,y,kR’l<I^u(x,y)u(x,y,kR,kj)dxdy  and 
w(x,y,kR,kj)  satisfies: 


1  ® ( x »y , kp , k » )  +  1 G(x-x  ,y* y  ,k_,k.) 


wlx.y.ko.k.)  =  e 


u(x' ,y' )w(x' ,y' ,kR,kj )dx'dy' . 

Multiplying  (20)  by  e"'il3(x»y»kR,kj)  and  employing  the  following 
symmetry  condition  on  the  Green's  function 


e-i^( x,y,kR.kj )G(x,y»kR,kj ) 
*  G(x,y,c0,kj) 


where  Cg  * 


—■k, ,  yields 
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whereupon  (10-11)  follow.  The  eigenfunction  ^  is  recovered  wit.n  rn* 
generalized  Cauchy  formula 

\  23 ; 


i  r rTf  (x >y 

u(x,y ,kR.kj )  =  1  ♦  7  fj - — -  dkR*dki 


noting  that  using  (10-11),  (23)  becomes  a  linear  integral  equation 
for  u.  The  potential  u(x,y)  is  recovered  by  taking  k-*-  in  (13)  or 
(14}  and  (23). 

For  the  K-P  the  evolution  of  the  data  obeys  ( v  *  4ik^  in  (9) 

H  ■  <8ik0)(6kk0  '  4k0  -  3k2)T 

o.k. 

where  kQ  «  kR  ♦  — ,  k  «  kR  ♦  ikp 

K 


Special  cases  include  o  =  oR-Hop 


(a)  KPjj  ;  o  3  -1:  oR  3  -l,Oj 


(24) 


|I  ■  81kR(3kf-k|)T 


(b)  KPj  ;  a  3  i 


aR-0-,  oj  3  1,  kj 


kI/aR 


(25) 


5T  ■  -81 ( kn-k, ) ( k|*2kRkj *4tf ) T 

These  formulae  allow  usin  principle  to  solve  the  Cauchy  problem 
for  K-P  and  in  particular  the  limit  (ii)  discussed  above  allows 
us  to  give  an  alternative  solution  for  KPj  via  5  and  not  via  a 
nonlocal  RH8VP. 

Similar  ideas  apply  to  higher  order  scalar  problems 

n 


(26) 


(iii) 

3v  .  3  v 
°3>  3xn 

where: 

v,  Uj  c  c  and  to 

o  £  • 

<jy  ax 

(iv) 

where: 

v,qitNxN, J=diag(  J 

,n-j 

,  ,3  v 


=  0 


Interested  readers  may  consult  refernce  4a,  and  review 
more  details. 

The  notion  of  1  extends  to  higher  dimensional  scattering  and  in¬ 
verse  scattering  problems.  However  as  we  shall  mention,  despite  the 
fact  that  the  inverse  scattering  problem  is  essentially  tractable 
there  does  not  appear  to  be  any  local  non  1  near  evolution  equations  in 
dimensions  greater  than  2*1  associated  with  multidimensional  gener¬ 
alizations  of  (iii)  or  (iv). 

Our  prototype  scattering  problem  will  be 

o  ♦  &v  +  u(x,y)v  •  0 
n  a2 

0*1  — *,  x  e  R  ,  y  r.  R  .  (27) 

t*l  3xJ 


Letting 


_  /  Jk*x  *  *  y/o 

v  «  u(x,y,k)e 


k  *  kR  ♦  ikj,  k  c  C 
k*x  «  l  k.x.,  o  •  c 

J  J  J 


°r  *  10r 


Then  there  exist  functions  u  bounded  for  all  x,  y  satisfying  u-1.  as 

|k.  |*  «,  j  *  l,...,n.  When  o.  I  0  u  turns  out  to  be  nonanalytic  in 
J  " 

each  of  the  variables  k,  i.e.  u  *  u(x,y,kR  ,...kR  ,  k.  . kj  )  and 

1  n  1  n 

satisfies  a  3  problem  linear  in  u,  in  each  of  the  variables  k.; 

J 

i.e.  we  shall  show  that  u  satisfies  an  equation  of  the  form, 


-ur-  *  T.(u);  j  *  I . n 

stj  J 


where  T.  is  an  appropriate  linear  integral  operator. 

The  basic  idea  in  order  to  derive  (28)  follows  a  similar 
format  to  the  two  dimensional  case  described  earlier.  From  the 
definition  of  u(x,y,k)  below  (27)  we  see  that  it  satisfies 


♦  iu  *  2ik-Vu  -  u(x,y)  =  0  . 


(29) 
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We  convert  to  an  integral  equation 

U  «  1  +  G(li  u  ) 

where  the  Green's  function  kernel  is  given  by 


i(x.£+yn) 

S - - - d£dy 


G(x,y,kp,k. )  ■  ''  *  2  „ 

R  1  (2*)  1  JJ  ioy-£  -2k*£ 


sanUL  _i_  [ 


.  Q  (-yoRU2+2  (kR+  ^-)fC)dC. 


Taking  the  3  derivative  of  (30) 


%  •  iP“’  * « «“  !*->• 


and  using 


3G/  .  i  \  1  i_  i  LiB(x,y,kD,kT ,5) 

ir(x,y,l‘vltI)  *  -  u^ri^1]6 


where 


• (C j"kRj )&(o(£ ) )d£ 


Slx.y.kp.kjC)  *  (x+2y  ^)-U-kR) 

R 


shows  that 


Ik-  *  -  — ^T  -M  T(ka.k.OUi-kRi)a(j(c)) 

3kj  (IT)  loR{  J  R  1  3  Rj 

•  w(x ,y *kR »k j , £ )d£ 
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where 


T(kRfkpf)  =  |  e  1-(x,y,kR,ki^u(x<y)  u'x,y,kRikj)dxdy  (37; 


and  w  satisfies 


/  ,  ,  j  i  S ( x ,y , k„ , k . £  )  . 

w(x,y,kR,kj,^  *  e  R  I  +  G(uw). 


*  1  fl 

fiiltiplying  (37)  by  e  and  using  the  symmetry  condition 


i S ( x  ,y ,  kR , k j ,5) 


G(x,y,kR,kj )  =  G(x ,y,C  *kj ) 


yields 


wfx.y.kp.kpf.  •  e‘i8<,,,,kR,kI,5)u(x,y,C,  k,) 


and  hence  (36)  gives 

5^  '  V“>  ’  -  T(kR.kI.e)Ccj-kRJ) 

•  6(cU))  ei6(  X,y,kR»lcI’C)u(x,y,c,kI)dc.  (4: 

We  see  that  "T  is  an  integral  operator  which  depends  on  a 
scalar  scattering  function  T  *  T(kR,kj,£)£  being  effectively 
(n-1)  integration  parameters  (due  to  the  delta  function  in  (41) 

in  the  nonlocal  operator  T.). 

J 

One  can  use  a  generalized  Cauchy  formula  such  as  (23) 
in  order  to  obtain  a  linear  integral  equation  to  reconstruct  u. 
However  due  to  the  redundancy  of  the  data  discussed  below,  we  find 
that  an  alternative  method  is  more  useful.  The  inverse  problem  is 
redundant,  i.e.  we  are  given  T(kR,kj,0  (3n-l  parameters)  and  we 
must  reconstruct  a  local  potential  u(x,y)  (n+1  parameters).  A 
serious  issue  is  how  to  characterize  admissible  inverse  data  T, 
i.e.  data  that  really  arises  from  a  local  potential  (small  generic 
changes  in  T(kR,kj,e)  cannot  be  expected  to  arise  from  a  local 
potential  u(x,y)).  Insight  into  this  question  is  obtained  by 
noting  that  T  must  satisfy  a  nonlinear  constraint,  one  which  is 
obtained  by  requiring  3^u/3k.9kj  s  3  u/sEjSE.  (i  t  j).  the  form 


-10- 


-11- 


trans forms  (42)  into: 


-~r  *  N< ;(T)(x,wn,w)  •  ja2,— n 
3x.  J  u 

J 


using  the  generalized  Cauchy  formula  (23)  we  have 

IjCTKx.w.Wq)  =  r(x*w,wQ)  -  |  -°2  dxRdXj 


*  u(w,w0) 


X  -  X 


46) 


(47) 


where 


x  ,  (X2.x3. — X j#— xn) 


u(wQ,w)  *  ||e"1^yVx'w^u(xty)dxdy  (48) 


We  have  used  the  fact  that  when  wQ  »  2kj . U-kR)/oR  and  w  -  c-kR 
are  kept  fixed,  T(x,w,Wg)  -u(w,Wg)  (The  Fourier  Transform  of  u(x,y)) 
for  large  Xj(w^fO);  this  is  the  analogue  of  the  Born  approximation. 

We  expect  that  for  suitably  "small"  u  (i.e.  no  homogeneous 
solutions  to  the  relevant  integral  quations)  if  I  is  independent 
of  x,j  and  decays  sufficiently  fast  for  |w|,|w0(*.,  then 
T(kR,kj,0  is  admissable.  Moreover  (47)  gives  a  formula  to 

reconstruct  the  potential  by  quadratures.  Limits  to  case 
3  s  i  and  reductions  to  stationary  potentials  u(x,y)  =  u(x) 
can  be  carried  out.  Oetails  can  be  found  in  Ref.  [5a, b].  It 
should  also  be  noted  that  in  recent  work  Nachman  and  lavine  [5c] 
have  extended  theaw  ideas  to  situations  where  there  are 
homogeneous  solutions  to  the  relevant  integral  equations. 

(42)  also  suggests  why  simple  local 


nonlinear  evolution  equations  have  not  been  associated  with  euuat’on 
(27).  Namely  in  theprevious  lower  dimensional  (2*1  and  1*1)  problems 
the  time  evolution  of  the  scattering  data  obeyed  a  particularly  simple 
equation,  (e.g.  ~  ■  w(kp,kj)T).  However  in  this  case  such  a  simple  flow 
will  not  be  maintained  -  due  to  the  nonlinear  constraint  (42). 


These  ideas  can  be  generalized  to  first  order  systems: 
(v)  ^  ♦  o  I  -  qv 


v.  q  c  € 


",H.  . jjj 

k  M  . 


with  many  similar  results  obtained  6a,b,c;  though  there  are  some 
important  differences  as  well:  see  ref.  [6c].  Again  the  scattering 
data  satisfies  a  nonlinear  constraint.  In  general,  there  is  no 
compatible  local  nonlinear  evolution  equation  associated  with  (v). 
However  when  certain  restrictions  are  put  on  Jj  then  the  constraint 
equation  becomes  linear  and  the  so-called  N  wave  interaction  equations 
are  compatible  with  the  system  (v).  Nachman  and  Ablowitz  [6a] 
showed  that  at  most,  the  system  would  be  3+1  dimensional,  and  Fokas[6b] 
showed  that  indeed  the  system  is  reducible  to  2+1  dimensions  by  a 
transformation  of  independent  variables  (characteristic  variables). 

In  [6c]  Fokas  studies  the  inverse  scattering  of  (v).  For  o  *  i  he 
finds  an  equation  similar  to  (42).  However  its  integrated 
form  shows  that  in  order  for  the  potential  to  be  reconstructed 
one  must  solve  a  reduced  system  of  equations  of  the  form  (v): 
i.e.  for  N  ■  2.  This  is  in  contrast  to  the  scalar  problem  where 
reconstruction  is  via  quadratures. 


Beals  and  Coifman  haven  an  alternative  but  similar  formula¬ 
tion  [7a, b]  for  multidimensional  scalar  problems. 

There  is  an  n-dimensional  problem  which  also  fits  within 
the  framework  of  1ST:  The  so-called  generalized  wave  and  generalized 
sine-Gordon  equation  (GWE  and  GSGE).  These  equations  arise  in  the 
context  of  differential  geometry  and  serve  to  extend  the  classical  re 
suits  of  Backlund  for  the  sine-Gordon  equation  to  n-dimensions  [8]. 
The  n-dimensional  Backlund  transformation  is  given  by: 
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dX  ♦  XA  a  =  A  -  XB, 


;  a?' 


where 


n 


v  ^X  . 

T —  <1*  i  - 

.  ,  <jx  .  1 

J’l  J  J 


dX  * 

Aij  =  6i(2)aijdXj- 
6 


and  a  «  { a,  J  e  F 

1  J 


nxn 


U  •  T77  srp  dxj  (50) 

Equations  (49-50)  reduce  to  the  Backlund  trans¬ 
formation  for  the  generalized  sine-Gordon  equation  (GSGE)  when 

S  i <2)  a  (22  ♦  (2«n  -  1)  )/2z.  (51) 

and  for  the  generalized  wave  equation  (GWE)  when 

S^z)  *  -(l-z2)/2z  =  x(z).  (52) 

The  compatibi 1 ity  condition  required  for  the  existence  of  solu¬ 
tions  to  these  Backlund  transformations  results  in  a  system  of  second- 

order  partial  differential  equations  for  an  orthogonal  nxn  matrix 
a  *  (a^j)  in  (49)  which  is  a  function  of  n  independent  variables 


a  -  a ( *  j .x^,. . . ,xfl), 


The  equation  has  the  form 


dX. 


_ 1_  ffn  3aik 

a  1  ka  1  j  3*k  9xj 


1  f  J. 


1  ,  j  ,  k  distinct, 


3aik 


a,  •  r»X 
ll  1 


i  *  v 


(53) 


where  <  *  1  for  the  GSGE  and  t  *  0  for  the  GWF. 
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We  observe  that  when  n  *  2  and 


1  (GSGE ) .  the  orthogona 


matrix  a  s  f  a . j  }  given  by 


cos  j  u 


sin  ^  u 


^-Sin  j  u 


cos  j  u 


for  the  function  u  *  u(x,t)  reduces  the  GSG£  to  the  classical  sine- 
Gordon  equation  ( <  *  -1), 


“tt  -  “  ‘  °- 


On  the  other  hand  when  n  ■  2  and  *  *  0,  then  with  64  )  the  GWE  reduces 
to  the  wave  equation  (55).  When  n  >  3  the  generalization  of  the 
wave  equations  discussed  here  is  nonlinear. 

The  BIcklund  transformations  (49)  described  above  are  In  fact 
matrix  Riccati  equations.  Linearizations  of  such  a  system  can  be 
performed  in  a  striaghtforward  manner.  Introducing  the  trans¬ 
formation 


X  «  UV_1. 


where  U,  V  and  n  x  n  matrix  functions  of  Xj. •••*„*  the  following  linear 


system  is  deduced: 


A\  / U 


with  the  components  of  A,  Q  given  by  (50).  Compatibility  ensures  that 
the  orthogonal  matrix  a  *  fa^l  satisfies  the  GSGE  with  (51)  and  GUE 
with  (52).  A1  ternati vely ,  if  we  call 


*  U'  . 


the  following  linear  system  of  2n  o.d.e.'s  are  obtained: 


=  AA.  4-  ♦  C  •*, 


iBUBraugaBKMSggiBgaaweueu^j^mjiniJWHnixuxiijniHtiix  its  ujhiexeu*  m  k»  mmgrwmn  w« 
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where  A.,  C  are  2n  x  2n  matrices  with  the  block  structure 
0  J 


(6C) 


Here  a-,  v  are  n  x  n  matrices  having  the  following  structure: 
J  J 

a'j  s  (r  ■  l)eidj  *  V 


aj  s  aej 


(61) 


where  e,  3  (e-)..  is  the  unit  matrix 

J  J  1 K  f 


{ej}  ik  m< 


l 


i  *  k  »  j , 


(62) 


and  in  component  form  y.  takes  the  form 

J 


0  otherwise, 
form 
3aai 


da. 


(Vkt "  (|-{kj>^;5if  4u  *  V631 

In  (61)  a  is  the  orthogonal  matrix  Fn  *  SO(n)  associated  with  the  GW£ 
when  6  »X  and  with  the  GSGE  when  6  *  ^(z  ♦  1/z),  X«  j(z  -  1/z),  and 
Yj  is  the  matrix  (63>:  Fn  -  Mn(F),  Yj  *  Yj  *  0.  Equations  (53)  arise 
as  the  compatibility  condition  associated  with  (58).  More  explicitly, 
for  the  GWE  the  scattering  problem  takes  the  form  [4)  *  4<(x,x)] 


sr'j  '“j*  *  v 


(64) 


wi  th 


(65) 


and  given  by  (60,63). 

For  tun  GSGC  the  scattering  problem  for  .  -.(»,:)  takes  the 

form 


A(Z) 


MRU  WWWKKWWWKWmrWWWW  11H«» 
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'■(*),  > ( z ) .  Cj  given  above. 


or  equivalently 


J 


■:  e  r  ■ 


where 

/O  ua  \ 

8j  5  I  JJ  ,  u  ■  diag< ♦  1 .  -1,  ....  -l).  (68) 

\a ju  0  / 

In  [8]  it  is  shown  how  these  linear  problems  may  be  vieweo  as  a  direct 
and  inverse  scattering  problem  for  the  GWE  and  GSGE.  Namely  the 
direct  and  inverse  problem  may  be  solved  for  matrix  potentials,  de¬ 
pending  on  the  orthogonal  matrix  a,  tending  to  the  identity  sufficient¬ 
ly  fast  in  certain  “generic"  directions.  It  should  be  noted  that 
solving  the  n-dimensional  GWE  and  GSGE  reduces  to  the  study  of  the 
scattering  and  inverse  scattering  associated  with  a  coupled  system  of 
n  one-dimensional  o.d.e.'s.  This  is  in  marked  contrast  to  other 
attempts  described  earlier  to  isolate  solvable  (local)  multidimensional 
nonlinear  evolution  equation  which  are  compatibility conditions  of  two 
Lax-type  operators,  e.g.. 


L  *  »  Ai(;  (69) 

4^  *  Mit  (70) 

where  L  is  a  partial  differential  operator  with  the  variable  t  enter¬ 
ing  only  parametrical ly .  Although  as  we  have  seen  nonlinear  evolution 
equations  in  three  independent  variables  can  be  associated  with  such 
Lax  pairs  (e.g.  the  K-P,  Davey-Stewartson,  three  wave  interaction 
equations,  etc.)  little  progress  via  this  route  has  been  made  in 
more  than  three  dimensions.  As  discussed  earlier  one  has  to  overcome  a 
serious  constraint  inherent  in  the  scattering/inverse  scattering 
theory  for  higher  dimensional  partial  differential  operators  in 
order  to  be  able  to  isolate  associated  solvable  nonl  inear  equations , 
i.e.  the  scattering  data  generally  satisfies  a  nonlinear  equation 
(eq.  (42)).  The  analysis  associated  with  the  GWE  and  GSGE  avoids 
these  difficulties  since  the  GWE  and  GSGE  problems  are  simply  a 
compatible  set  of  nonlinear  one-dimensional  o.d.e.'s. 
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The  results  in  ref.  [8]  demonstrate  that  the  initial  value  problem 
is  posed  with  given  data  along  lines  and  not  on  (n-l)  dimensional 
manifolds. 

Similar  ideas  apply  to  certain  n-dimensional  extensions 
of  the  so-called  anti-self-dual  Yang-Mills  equations  (SDYM) 

[9].  In  two  complex  variables  the  self-dual  Yang  Mills  equations 
take  the  form  (see  [10]} 


4-  (G-1  i2_ }  *  JL-  (a* 

3xr  3XJ  ax2 


1  !£L)  . 


)  *  o. 


where  &  is  a  positive  matrix  valued  function  of 

Alternatively  SOYM  takes  the  form 

3A.  3A, 

_L  +  ,  o 

3xi  3*2 

3A,  »A, 


4  -  4  *  •  0. 


where 


A  .  »  -n 

J 


-1  an 


The  SDYM  may  be  obtained  via  the  compatibility  condition 
of  the  following  linear  system 


3m 

-  z  -7T-  3  A.  m 
®X-  A 
2 


a*2  M  ‘ 

multidimensional  extensions  may  be  obtained, 
consider  the  linear  system 


0^m( x ,z )  »  Aj(x)m(x,z),  j 


*  —  ♦ «,  4 — 

"j  J  3*ri 


For  example. 


1 *  n 


Vl  *  X1  •  sj  '  (-UJ. 


iv 


h. 


k- 

V;* 

1  v 

iv< 


to 

KS 


$ 

$ 

£ 

I 

$ 

k 

w 

& 


ss 

I 


ft 

i 
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Compatibility  (commutativity)  implies: 


°iAj  -  DiA.  ♦  tvv  *  * 


3A.  3A. 

— - rJ-  ♦  [A,, A.]  »  0 

3xi  L  1* 


A  potential  ft  may  be  introduced  as  before: 


»  .  n*1  2S_ 

i  »«j 


to  obtain 


.  3  /„-!  3ft  \  3  3Q  \  .  n 

sjl^(n  *  *1  T^>*° 


Clearly  when  n»2  this  system  reduces  to  the  classical 
SOYM  equation. 

Solutions  to  these  equations  may  be  constructed  via 
the  3  method.  Define 


D-j  ■  l\  *  zlj 


d  I  J  *  e  _ 

3x.  ’  2  j  , 

J  J  +  l 


We  shall  show  that  the  3  integral  equation 


m( x  ,z )  *  I  ♦ 


mV)(*, 


dc-d; 


»  ri  i'i  i'i  i't  t 


i  .*1  ,'|  *'| 


|l|  .fi  ,|i  . 


satisfies  (76).  Operating  on  (84)  with  0^  yields. 


.  f  (L^m)V  ♦  md-jv) 

)Jm  *  t*7  — - - — dc-dc  ♦  J 

2  2,1  J  C-2 


where 


zLi(mV) 


in  I  ;-z 


d;*dc 


a  ‘  jk  \  L^mV)dl;‘dC 


,  (•  cU(mV) 

Ar  j  -tt- 


Putting  (85).  (86)  together  gives 


.  .  ,  f  (Djm)V  ♦  m(0jv) 

"J-  ■  V  J  8  *•« 


where 


Aj(«) 


U  lj(mV)dc-dc  *  -  s.  -rj. —  f(mV)dc*d5  .  (88) 

ITlJZ  J  +  l  j' 


j  «j  +  l  1 


We  shall  require  V(*,z)  to  satisfy 


D^V  «  0 


in  which  case  using  (84)  in  (87)  by  writing 


Aj  ’  V"  -  ITT  f?idc'd*) 


ikV 


[v 


Fi  ' 
:< 


* 


»» 


k: 


K: 


© 

y 

: 


h?: 


I 


(89)  v 


1 

Ml 


£ 


(90)  S' 


i 


tv: 
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we  find 


.  .  }r  - 

(0j"  -  Ajm)  ’  27T  J  - Z 


((0 ;m)  -  A.(x)m)V 


“d^ -  dr . 


!Tf  i 


For  y  suitably  chosen  (84)  has  a  unique  solution  in  which  case 


O^m  -  A  .m  •  0 

•  w 


Thus  A.  «  A.  and  solutions  of  the  extended  SDYM  are  obtained. 

J  J 


The  condition  (89)  is  satisfied  if  we  take  V(x,z)  *  V(u(*)»z). 
with  Uj(x)  »  zxj  ♦  +  i  and  V  holomorphic  in  the  Uj.  Then 


°iv  ■  <3ib  *  »<v—vl) 


*  SjSjHZ>>  =  ° 

by  virtue  of  s^  «  (-)J.  In  ref.  [9]  other  examples  of 
multidimensional  extensions  of  SDYM  and  a  rigorous  derivation 
of  the  foregoing  is  given. 
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